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SUMMARY 


c 




By an. extension of Eward's "diaphragm" concept outside the wing 
tip, the present paper presents two approximate methods for calculating 
the aerodynamic behavior of harmonically oscillating, sweptback finite 
■v/lngs with both supersonic leading and trailing edges. Both methods 
are based upon the fact that the contribution of the so-called "diaphragm" 
on the potential at a fixed point on the wing is approximately canceled 
by the contribution of a portion of the wing itself, if the wing under- 
goes steady harmonic oscillation at a low value of the frequency param- 


eter 


cix: 

ap^ 


(where cu is angular velocity per secondj c is chord; a 


is velocity of sound; and p = - 1, M being Mach number) . It is 

further necessary that the point be influenced. by the diaphragm from 
only one wing tip. 


The first method expresses the integration in terms of Fresnal 
integrals. The final expressions for the presstire coefficient are in 
the form of integrals which may be readily evaluated, • using numerical 
integration with available mathematical tableg. Some calculations are 
shown, and the results check well with the results of the second method. 
It is expected that this method should give accurate results in the 

engineering sense if the 'frequency parameter — 5 - is of the order one- 

aP'^ 

half or less. Even if the frequency parameter exceeds this value, the 
pressure increments calculated asymptoticEjJy approach the correct values 
at the inner edges of the tip Mach cone. 


As a byproduct of the above approach, the expressions for the down- 
wash of the wing wake, influenced by only the two tip diaphragms, are 
given in terms of, Fresna l integrals. 


The second method is similar to the first one except that a series 
expansion of the soizrce- strength function is used. When the frequency 

'ok ^ 1'' 


parameter is small 


< 


ap' 


the few terms used give good engineering 


'T- '' 

fiKl 
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accuracy. In the present case, five terms are considered. Of course, 
the numher of terms must increase as the frequency parameter increases. 
This method gives closed expressions for the pressure coefficients. In 
applying the method, the calculations are lengthy, hut no great cnmpnt lng 
skill is required. For the benefit of design engineers, some expressions 
of the average pressiure coefficient over the entire rectangular wing tip 
influenced by three-dimensional flow are given in a closed form, so that 
they can be used to estimate the tip effect in preliminary design. 

A number of graphs are presented to show the pressure distribution 
of the wings at different I4ach numbers, sweepback angles, and frequency 
parameters. Coinparlsons of the pressure coefficients calculated from 
the present methods and the exact two-dimensional solutions based on 
linearized theory are made. 


HJTRODUCTION 


The aerodynamic behavior of a harmonically oscillating finite wing 
in a supersonic flow has been recently studied by a number of investi- 
gators, particularly Garrick and Rubinow (reference 1), Eward (refer- 
ence 2), and Miles (reference 3 )- From the viewpoint of mathematical 
physics, it is one phase of the mechanics of wave motion investigated 
by a nmber of great physicists and applied mathematicians such as 
Huyghens (reference 4), Maxwell, Kirchhoff, Sommerfeld, Hadamard, 
Volterra, and others over a number of years. The classical treatment, 
however, is principally along the line of Cauchy's initial -value problem, 
while the present aerodynamic problem is a boundary- value problem. The 
classical resiilts are of some help in the formulation of the present 
problem but not of much help in the solution. 

With the advent of griided missiles and aircraft operating in the 
supersonic range, trouble has -been experienced with flutter of the wing 
or control systems. It was such troubles that spurred the author to 
obtain an understanding of the aerodynamic behavior of a harmonically 
oscillating finite wing at supersonic speed. The two-dimensional case 
of an arbitrary unsteady motion has been solved in the last few years 
(references 5^ 6, 7^ and 8) . But in the three-dimensional case i*p to 
now very few useful results have been achieved. 

The present paper treats the case of a finite wing with the leading 
edge ahead of the Mach line. The wing may be swept back or forward. 
Another condition has to be Imposed, that the Mach cones at the nose and 
the tip never Intersect within the wing area. With slight modification, 
the method can be extended to more general wing plan forme. The wing 
is s^q)posed to undergo haraonlc oscillations of small amplitude during 
flight at supersonic speed. 
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For the wing plan form as specified above, the three-dimensional 
flow fields are limited to three portions of the wing area: the area 
within the nose Mach cone, and the two areas within the tip Mach cones. 
In the remaining wing ai^eas, the aerodynamic behavior can be calculated 
with the two-dimensional theory of swept wings with supersonic leading 
and trailing edges under harmonic oscillation (reference 6, e.g.). 
Therefore, the present investigation is limited to the tip area and nose 
area of the wing. 

The author should express his appreciation to Miss V. O’Brien, 

Miss P. darken, and Mr. B.. T. Chu for their assistance in carrying out 
the project. The present project was conducted at The Johns Hopkins 
University under the sponsorship and with the financial assistance of 
the National Advisory Committee for Aeronautics. 


LBJEAEIZED DIFFERENTIAL EQUATION AND 
SOLUTIONS. OF SOURCE INTEGRATION 
Genereil Analysis 


The linearized differential equation of unsteady compressible flow 
is well-known (references 1 and 5) 


3x2 5y2 a2\3-t 3x/ 

or, briefly. 



(la) 


(lb) ■■ 


where U • is the free-stream velocity and a is the sound speed of the 
free-stream condition. (Symbols used in_thls paper are defined in 
appendix A.) The velocity vector q = Ui + q' = Ui + grad 0 and 
is the Laplacian operator. The elementary solution of the potential at 
a point (x,y,z) at time t due to a source of time- dependent strength A 
is 




~ ^l) A(|jT],^,t - Tg) 


r 


(2a) 
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■wiiere the soitrce is located at and 


r = - i)^ - 3^[(y - + (z - 


( 3 ) 


_ (x - £)M r 
ap2 a 


T = 


a3‘ 



W 


where 

- 1 


In the case of harmonically oscillating source strength the source- 
strength function can he written as a product of two functions, one 
containing the source position and the other showing the time dependence, 
or, explicitly. 


A(l,Ti,^,t - t) = A(|,Ti,^)f(t - t) 

Equation (2a) then hecomes; 

5,TJ,0 = jf^t - + f^t - (2b) 

In the whole domain surrounding the wing in supersonic flight, there 
is one particular streamline surface which is of interest to the problem. 
This streamline surface is formed by the Btreamlines, which coming from 
upstream either strike the wing leading edge or pass throu^ the fore- 
most' Mach line of the tip Mach cones. It follows the contour of the 
wing (both upper and lower srirfaces), reunites at the trailing edge, and 
then extends to infinity at downstream. 
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Figure 1 shows the streamline surface which Is associated with an 
oscillating wing In a supersonic flow- at time t. 

The whole streamline surface can he divided Into eight zones as 
shown. Ahead of the supersonic leading edge and outside the two tip 
Mach cones, the flow Is not Influenced by any of the disturbances caused 
by the *wlng. This Is called zone 1 as shown. The portion of the stream- 
line surface behind the supersonic leading edge and between the tip Mach 
cones Is called zone 2. In this zone the potential at any point on the 
top wing surface (or bottom) Is never Influenced by the flow on the 
opposite side. Therefore, the source solution of equation (2b) can be 
applied without restriction. If the source-strength function A Is 
known. 


Zone 3 represents the portion of the wing Influenced by the right 
tip Mach cone and zone by the left tip Mach cone. In these two zones, 
the flow on one wing surface Is Influenced by the flow on the opposite 
side. The sli^ple source solution of equation (2b) cannot be applied. 

The present paper gives the approximate solution for these two zones 
for the case where the two tip Mach cones' meet on or outside the trailing 
edge. The present method with slight modification can also be extended 
to the case that the tip Mach cones Intersect Inside the wing area (not 
shown In fig. l), but no diaphragm within the forward Mach cone of the 
point under consideration Is Influenced by the dl^hra^ on the other 
side . 


Zone 5 Is the portion of the streamline sijrface outside the wing 
surface Influenced only by the right tip Mach cone; similarly zone 6 Is 
Influenced only by the left tip Mach cone. Zone 7 is the portion of the 
streamline surface after the trailing edge Influenced by neither Mach 
cone. Zone 8 Is the portion of the streamline surface Influenced by both 
Mach cones. Of course. If the wing moves slightly up and down with time, 
the streamline surface and consequently the diaphragm must also move up 
and down with It. With the concept of Eward (reference 9)^ any finite 
wing in supersonic flow can be considered to be attached over Its stream- 
line surface outside the wing (such as zones 5^ 6, and 8) to a thin 
weightless rubber diaphragm, of which the pressure on both sides at a 
point must be balanced and throiigh which no fluid can leak. 

The diaphragm Is imagined to extend to the outer edge of the tip 
Mach cone, and to coincide with the shape of the streamline surface. 

Thixs, the rubbfer diaphragm can be considered as an extension of the wing. 
The flow on the upper side of the wing with Its extension can be con- 
sidered as Independent of the flow on the lower side of the wing 
vice versa. The potential at any point on one side of the wing can be 
evaluated from the soiarqe distribution on the same side of the wing and 
the diaphragm within the forward Mach cone of that point. In the steady 
case It iias been shown by Puckett that the source strength at a point 
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of the wing is directly proportional to the slope at the point and inde- 
pendent of the slope of its neighhbring points. Eward has shown that 
the solutions obtained using the diaphragm concept are identical to 
those of conical flow (reference 9) • 

For steady flows, the principal advantage of Eward's method is 
that the slope of a point on the diaphragm is determined solely by the 
slopes of the wing surface along the characteristic or Mach line through 
the point, as long as the point on the diaphragm is ahead of the Mach 
cone from the other wing tip. In the case of unsteady flow, it can be 
shown that the above statement is approximately true if the wing slopes 
vary very slowly with time.^ Using the diaphragm concept, the disturbance 
potential at any point can be considered as the effect of sources of the 
wing system on the same side. (Of course, the slope of the imaginary 
diaphragm depends on the flows over both wing surfaces.) It becomes 
apparent that the potential can be evaluated by the well-known source- 
integration method. Let us consider a point on the diaphragm where the 
potential contributed by sources on either side of the wing must be the 
same. It has been shown by Eward that the sotirce strength is propor- 
tional to the vertical component of the perturbation velocity, or, more 
directly, is proportional to the local slope of the wing surface or the 
diaphragm even if the flow is not steady. Mathematically, this is shown 
w Uoc. • 

as .A = — o = — o where w is the vertical con 5 )onent of the pertur- 
bat ion velocity and a denotes the slope of the streamline surface. 

The following equations can be written to denote: 

The elope at the upper wing surface at time t - by 

= Oiji^l^'n,+0,t - 

The slope of the tipper wing stirface at time t - T 2 by 

af,2 = " '’’ 2 ^ 

The slope of the lower wing stirface at time t - by 

%1 = - Tp) 

The slope of the lower wing sttrface at time t - T 2 by 

%2 = " ^ 2 ) 

The slope on the diaphragm at time t - by 

= ^(l^n,0,t - 

The slope on the diaphragm at time t - T 2 by 

^ - ^2) 
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Then the potential can be written at any point on the top or bottom 
streamline sinrface as: 


^(x,y,40,t) 






+ X2) 

r 


-0;,t) 



(‘^,1 ° B , 2 ) 

r 




( 6 ) 



Or^ in an oblique coordinate system whose axes lie parallel to the Mach 
lines and whose origin is placed at the point of tangency of the Mach 
line to the leading edge of the wing tip (the junction of the supersonic 
and subsonic leading edge) (see fig. 2 ), these equations can be written 
as: 




where (u.v) is the transformed coordinate of the point originally at 
(x y) and (u',v') is the transformed coordinate of the soiurce eluent 
oJiglnSy at^M) in the x,y system. The old ^d new coor^te 
systems are related hy the following transformation eq^IationB. 


u = ■^ (x - Py) 


T = (x + Py) 


X = £ (v + n) 


y = i (v - u) 


u' = (I - PTl) 


V' = ^ U + Pn) 


I' = I (v’ + u') 


Tl =1 (v’ -_U’) 


d| dn = % du' dv' 





N9 


and tlie quantity r and the time intenralB and Tp in the new coordinate Byatem are_, 

respectively: 


f - \J(u - u')(v - v') 


""i “ ^ ^ - ilid 

° i + u - u') + ^ \/(u - u')(v - v') 




(8b) 


If the supereonic and subsonic leadiiifi edges are denoted by the equation (see figs. 2, 3, and 4) 

V = V 2 (u) or u => U]_(y) 

V = V 2 (u) or u = ti 2 (v) 

then the potential at any point Vv) on the wing can be written as: 


(9) 


u 




du' 




^ - "^l) + - T2) , , 

r .w ' g ■ 


fuo(vv) 

du' 

- u'j 


'->v^\u- j 


(Jjl 

u',v',t - 

' ^1) 

1 + CTj(u’,v',t - T 2 ) 

1 

s- 

, v>; 

(1/2 


1 dv' + 


^T2(u') 


x(u',v',t - Tl) + x(u'jV',t - Tg) 


(10) 


VO 


19ns Ml VOlfH 



(a similar expression can "be obtained for ^.) Now it will be shown that in the case of har- o 
monic oscillation the second integral 'is approximately zero under certain conditions. 

For steady harmonic oscillation of a thin flat wing, assume: ^ 

-Og “ aip = a a cLe-*^ (lla) 

where a is everywhere the same on the wing. 

Further, assume: 

\ = X(u,v)ei^ (11b) 

where X(u,v) is the amplitud.e of oscillation, of points on the diaphragm and is a function of 
the point (u,v). The term X(u,v) may be complex, including the time lag of the oscillation, 
because the time lag must be a function of location only. 


With a and \ given by equations (lla) and (lib), respectively, the second integral 
(along constant u) of eq'uation (lO) can be written as: 



Sm VOVH 





When CD = 0, that is, the steady case, this integral reduces to Ervard's result in refer- 
ence 2 and has been shovn to be exactly zero. For co ^ 0, the integral ■will again be zero if 

COB (u„ - u')^/^(tv - 

fipH 

consideration. 


= 1 exactly. This will be clear from the following 


Eward in reference 2, equation {2k), shoved that for any point 
(see fig. 2) the following relation exists along the characteristics 


^Up, v-Q^ on the diaphragm 
line u' = Up! 



^ o’B(^Sv',t - T]^) + Ug(u',y',t - Tg) - a^(u’,v’,t - t^) - aj(u',v*,t - xg) 

/ ' ~i ' . . T /O ‘ ^ 


J. 


Vi(u’) 


. . T /O 

2(Vd. - V'p - 


<iv' (13)' 


This means that along the characteristic line u = Constant the total effect of the - diaphragm 
source strength is equal to the total effect of source strength on both wing surfaces along the 
same characteristics, and that it is Independent of the source strength of the diaphragm or the 
wing surface outside the above characteristics. Now, along this element u’ = Up, 


^Members of the HACA staff recently expressed the opinion that "Since and Tg con- 

tain the variables u', a second-order integral eqiiation is generally requljred rather than the 
simplified first-order Integral equation given. Equation is satisfactory if t may be 

represented as varying linearly with time in the interval (^i' “’■’' 2 ) • ” ^ present analysis, 

the use of equation (13) is Justified. 


JU VOViJ; 
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= ^ (’'D - ^') = 


SubBtitutlug equations (lla), (lib), and the above value of the t’b 
into equation ( 13 )^ there is obtainfed: 


4-^) 

X(u',v»)e ' ^ ^ 

(vd 


4g(u') 


iCD t 


f2(u' ) 


- 


dv' = - 


oe 


(-¥) 


^Vt(u') 


(VB - vf 


dv' (ll<-) 


Nowif in equation (12) cos ('‘^w “ ~ "v y • 


r>\l/2 


= 1 exactly. 


the Integral at the right-hand side of the equation will reduce to the 
form: 


'Tg(u') ^ 




w 


-ln> 


iTi(u' ) 


K 


v')V 2 


dv* + 


Vu') 


X(u* >v' )e ^ 

(■% - v')V 2 


N-v') 


dv' 


(15) 


which according to equation (l4) is zero, beca^^Be here v^ = Vp 
In general, cos ~ -.1^1/2/^ _ „t\l/2 


,"w 


- - ■^T'" 


cnr 


= cos — is less 
a 


than unity. (Note that r = -^ ^(x -_i)^ - P^(y - is the hyper- 


bolic distance, not the geometrical distance.) But for smal 1 values of 

S\q)pose that 

Then 

a — c; 

2 


car 

a 


CUT CUT 

— , the variation of cos — from unity is very small 

6L Q. 

the following condition is imposed: — = o* 

€L ^ 

1 '2 COB ^ ^ 1 - i (^) = 0.875' In both integrals the factor cos 

will reduce the value of each Integrand a small percentage in view of 
the monotonicaH ly 'decreasing nature of cos — for small values of — . 

S, €L 

Of course, the net contribution of both integrals cannot cancel out 

CDI' 

exactly as if cos — = If which is located on the left tip Mach line. 



ft 


a 


For a first appro xim ation , however, there can he written that the contribution of the 
diaphragm Sj, (fig. 2) and the portion of wing area cancel each other; otherwise, the 

aolutlon 1b too difficult to obtain, particularly with aweptback wings. Equation (lO) will 
then become; 





du' 


(% - 




0|pl 

-p 

> 

- Tj) + 

[u',V’,t - Tg) 

1 


,1/2 


dv’ (l6) 


If the above approximation Is analyzed more carefully. It will be found that the average value , 

(ret* 

of cos — on the diaphragm 1b slightly larger than that on the wing. This can be shown In 
the following way; On the diaphragm V 2 ^ v' ^ v^, and on the wing ^ V ^ Vg (fig. 2). 

For the purpose of Illustration, take the average value of v' =■ on the diaphragm and 

Vi + Vg ^ 

= on the wing. Then it is clear that 



^ ^ o'ther words, the contribution of the diaphragm is slightly more than the 'oppo- 
site contribution of the 'vdng area to the tip corner of the trailing edge, iJartictaarly 
on the Bide edge of the tip -where the contribution due to the wing area Is zero (see 

equation (l6) and fig. 3). Therefore, the pressure at 'the side edge of the tip cannot be zero 
as predicted with -the present approximation. But as the point considered approaches the left 


1^3 HI VDVH 



I 


I 


characteristic of the right Mach cone (l.e., r = O) the preseat approximation approaches the 
exact solution as a limit. This la faverehle to the present approximation for two reasons. 
First_, the pressure Increment In the nelghhorhood of the side edge of the tip is e::qjected to he 
negligibly amall^ particularly at a low frequency parameter ■ Its contribution to the aero- 
dynamic behavior of the wing specified is negligible. The present approximation asyn^itotlcally 
approaches the exact value near the left tip Mach line where the magnitude of pressure incre- 
ment is much higher and becomes more lD 5 )crrtant to the aerodynamic behavior of the wing than 
that at the side edge of the tip. Besides, the major portion of the wing area is calculated by 
the tvo-dlmenslonal theory as shown in references 5 and 6 and total contribution of the tip is 
usually not too important. Therefore, the approximation Is of good engineering accuracy. 


Calculation of Dovnwash Angle X.(^u,Vjj,tj for Harmonically Oscillating Wing 

In- reference 2 (equation (26)), Eward has shown that the downwash angle in zone 5 or 6 
(fig. l) is given by: 



( 17 )^ 


An approximate method of evaluating the above integral is presented below. 

If equation (lla) is substituted into equation (17)^ the following equation is obtained: 



%ee remarks about equation (l3). 
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Obviously, 


' X — > 00 as Vj) — > V2 
X — » 0 as Vj) — 00 


For the case Vp > Vq (fig- 2 ), V2 - v^^ can be divided into n equal 
parts and the 1th division can be denoted by Vjj^' (where Vj^’ = Vj^ 
for 1 = 0 and v^’ ~ '^2 1 = n) . Then, 



where the division is so fine that V2 - v’ = V2 - v^’ can be taken out 

of the integration sign for the first approximation. Incidentally, the 
integral left can be transformed to the Fresnal integral (reference 10 ). 
Equation (I9) then becomes: 




l6 
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Some of the basic properties of the Fresnal Integral are given in 
equations (29) • 

APPROXIMATION THEORI OF VELOCITY POTENTIAL AND PRESSURE COEFFICIENT 
AT A POINT ON A HARMONICALLY OSCILLATING SWEPTBACK' FINITE 
WING WITH SUPERSONIC LEADING EDGE^ - FIRST METHOD 

Asstmie the vhole wing flapping up and down according to 
Vip(t) = 06'^^ where a is the an5>lttude of oscillation. In fact, 

_ Maximum descending rate of wing hjuax 

a = 

Free-stream velocity U 

The frequency of oscillation per second is — . The velocity potential 

2rt 

of such oscillation may be considered as due to harmonic sources and is 

^The author should express his appreciation to Dr. J. C. Eweird 
for a discussion on this approximation method in May 19^. 






given "by equation (6) . It lias also been shown that for the harmonically oscillating wing, 
equation (6) can be further almplified. to equation (l6) as follows: 


Doe 


icut 




-ioyn -ICDTo 
e + e 


dT\ 




U 

' 2 m 


P 


% 




du* 

r"" 0,, 

fu’,v’,t ■ 
> 

- Tt_] 

1 + ( 

Jrrifu’,v’,t - To) 


|l /2 i 

f 

1 

(-W- 

. V) 

, 1/2 


where 


qip(u',v',t - Tjl) = oe 


Im 


(t-Tl) _ 


Ifl) 




= oe 


o^(n',v^,t - Tg) = 


_ 103(t-T2) _ icD t -^(vv-v'+uv-u’) -^(uw-u')V2jy^_Tr^l/2 

oce = ote 


i^ns tSL VDVK 
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Since in most practical^ cases straight leading edges are under consider- 
ation, put (cf. equation (9) fig- 3): 


V = V2(u) = kj^u or u = U]^(v) = v/k^ 

V = VgCu) *= k2U. or u = u^Cv) = v/k2 


as the equation of the supersonic wing edge outside the tip Mach cone 
and the subsonic leading edge within the Mach cone. Then, equation (22a) 
cein be written as; 


Ji(u„,v„,-K),t) - - 


Ude 


2Mtt 



du' 


("w - 


372 


(V. 


dv' 


(^v - 


:^[(v„-v-+n„-u' ) +|(«„-u' )^^(Vv-V ) V2] 

372 ® 


k.u 


where it is understood that e 
stands for the sum of e 

^ [( Vv-T' +u^-u' ) + |(uw-u' ) V2(v^-v‘ )l/2j 

and e 

Introducing the following new variables Z and Z instead of v', 

2 - i 

so that 




Then eqiiatlon (22a) can he written ae: 


Uoe"" /2 cq\ 


du'e 




PZ:. 






iZi ?2 

" 2 ^ 

dZe + 






where 


Uae^/2cD ^“V2 




M3t W 


»)- 


du'e 


aM^ 


J 


^wA2 


(Sf - 


-C(Zi) - c(Zi) + l[s(Zi) + S(Zi)]l (24) 


2tD\i/2r. 


1 ^ lapn/ “ ^1^') 


Zi p 

\apji/ l: 




correaponding to v' 
cortreaponding to v' 




= V. 


_ ,^\l/2 
^1 - U 


- kiu')V2 + i _ u')V^ corresponding to v’ 


w 


= k^u' 


2, , . „.)V2 


M tort/* 


corresponding to v’ = v, 


w 


(25a) 

(25b) 

( 25 c) 

(25d) 


>10 


IS^lZ Hj; VOVH 
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Wow introduce another new variable Y instead of u’ : 



80 that 


1 /2cdP \ 1/ 2 du» 

^ InaM^/ (ti„ - u')^/^ 



2 to . 
afiit 


^1^) + 



Y 


Z 


2 




2m 

aPn 


(^w - + 


kiM‘= 


1/2 


3 


(26) 


(27a) 


(27b) 


(27c) 



(27d) 


If 


Yl = 


/ m W/g, 

\naM^/ 


^ - 
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then equation (24) becomes: 


0(^wVw+O,t) 


Uoa 

(a 



0 


dYe 


irt„2 

2 


-C(Zi) + IS(Zi) 
-C(Zi) + iS(Zi) 


(28) 


At this stage more Fresnal integrals can be introduced by integrating 
by parts. For this purpose, the following properties of a Fresnal 
integral may be useful; 


C(e) - iS(0 = 




1 

2 , 


e d.-=| 

0 ^0 

lrt€£ 

2 

H_i/2^2)(^.) 


Ojte^ 

2 




da' 


(29a) 


(a. = I 


S(e') d6' = €S(e) + i cos £ - i 

n 2 n 


(29b) 


C(e) 

S(e) 
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■where j( ) is the Bessel function and ) is the Hankel function 

of the second kind and various subscripts indicate the order of these 
functions. Por large "values of e. 


“ I TO I °(^) (29f) 

(29g) 

s(±«) = C(±«>) = ±“ (29h) 

S(e) = -S(-e),, C(e) = -C(-c) (291) 


For details j see reference 10, pages ^kk to 546 and 744 to 745, and also 
reference 11. 

However, it is difficult, if not in^jossible, to reduce equation (28) 
to a closed form. Besides, there is usually more interest in the pres- 
sure coefficient which is generally kno'wn as: 


Cp(x,y,t) 


A M _ 2 M 

u2 ht U 8x 


(30) 


rather than the potential itself. In the u,.v coordinate system, the 
pressxire coefficient is gl-ven by: 


Cp(x,y,40,t) 


tj 2 8t u\8uw 8vyr 8x/ 


_ _ 2 _ ^ M ^ 80 V 

u2 8t UP\8u^ ^ 8v^/ 


(31) 


' Now^ 


8t 


= hD0(u^,v^,+O,t) 


(32) 
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Differentiating equation (28) with respect to and v, 




Irt V 2 
— ;r 


e ^ ^ 


Y]_ 


-C(Z3) - c(z^) + i|s(Z3) + s(Z3^ 


Ukl 5 e^ 


.i £,2 

e ^ dYl^r- + 


PZ^^ + Y PZ^ - Y 


-M. = __L UP5e^^ e 


- — Yt ^ 

-.icut 2 1 r- 


Yj_ 


-°(^ 3 ) - <=(%) + ^(^3) * ^(^3) 


Uae 


Ixiyb 


nil 

e dY F + 


rt 


PZj + Y PZ2 - Y 


where 


Z3 - 


2a> , \ kiM^ P 

^ (^yr - >^1'V) ■+ il^ 


1/2 


Yl 

^ 


Z3 = 


2 o) p 

(Vv - kl-w) - ^ Yl" 


-11/2 


(33) 


( 31 ^) 


(35a) 


(35b) 
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The pressure coefficient is then given by; 


C-D(x,y,40,t) 


' [''("3) - <^("3) * iS(Z3) 


^Imt/ 1 A - M 1 ^ " 2 


iS(Z 3 )] 


PY- 


1 i« v2 ^ 

“2 

dY e <; 


-f - ‘=(^ 1 ) + * 




is r/ 2 is 7 2> 


PZ, + Y PK - Y 


( 36 ) 


The above equation is used to calciaate Cp along the trailing 

edge of a rectangxilar wing tip which is harmonically oscillating. (See 
fig. 5-) In this particular Instance, M = \[2f the frequency parameter 


OK 


- 5 ' = 0.2, and the chord c is teiken as - un i ty. The curve is nlotted in 

the form — against y (positive y measured outward along the 

— icut ' 

ae 

span) so that for any specified amplitude of oscillation a the pres- 
sure distribution at ai^r instant t can be estimated from a single 

Cp 

curve. As is a congplex quantity, the real and imaginary parts, 

denoted by j ^ respectively, are plotted separately. 

If the harmonic oscillation is a cos ost, that is, the real part ' of 
oe^^, then Cp can be calculated from; 




a 






cos Ojt - 


;^=ioJt 


iOe 


sin cut 


(37a) 



t t - - r 


Similarly, if the oacillatlon la a sia rot, that 1 b , the imaginary part of oe-^^, it 1 b given 



uy : 



a 



aln ait + 



cos rot 


(37b) 


It 1 b Interesting to note that when y = 1 these values conform with the two-dimensional 
values given In reference 6. 


As soon as Cp everywhere over the wing, tip 1 b known, it Is simple to obtain Cjj, 

and Cy graphically. ■ 

As far as the calculation of the center aection 1 b concerned. It is easy to show that 
(see fig. 4) 




r 


TiW 


Pv, 


airt 


I 


w. 


dv’ 


(y - 

Ukg'u' ^ 


Iroft-Ti) lillft-To) 

Q— ' + e ' ■ 


TJg 

2Mrt 


PO 


du' 




v„A,' ('V - 


‘dv' 


uki 'u 


•nW ('W 


v„ - V 


') 


1/2 


^lro(t-Ti) ^ ghn(t-T2) 


(38) 


ro 


l9^3 HI ViJVH 



where 


ro 

os 


1 

■^1 “ ap 


^w - ^ - I ("V - 


■"2 ” - ^' + ^ - ^' + i (^ - 


Introducing Z and Z as "before (cf. equations (23) )^ and then following the procedure 
of deriving equation (24-), equation (2o) can "be written as: 


rf(u.v.+0.t1 

I % -- ^ ■ f f r 


use 


lent 


! 2od 2 


I ^ aM^ I- 


nM \aPny 


du' e ' 


J, 


(“» - 

^Uv-u') 


-n/z-A _ n/’7„'^ 4 - iR/'7.. 4 . ih/7.. 

“\“JL j ''\"l. / ■ } ■ — \-± / 


p° ■ -iu:p 

^^wAi 


where, replacing k3^ "by k2*. 


Zi’ = Z]^ in equation (25a) 
Zj_’ = ^ in equation (25c) 
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and, replacing by k^'. 


Z]^" = Z^ in equation (25a), 
Zj" = Zj^ in equation (25c) 


Similarly Y can be defined as in equation (26) and by tbe same 
procedure Cp can be determined. 

APPROXIMATE ARAIYTIC SOLUTIONS OF PRESSURE ON AN 
OSCILLATING FINITE WING - SECOND METHOD 


In the preceding section, an approximate integral solution for the 
oscillating finite wing is given, and an example of calculation is shown, 
but the full evaluation of the integral for a large number of cases 
would take more time than was available. However, for very low fre- 
quencies, this problem can be solved asymptotically by using a few terms 
of a series expansion of w(u,v,t). In the following section the approx- 
imate a n al y tic expressions will be shown. “ ' 


Analysis 


From equation (2^a),' 


05,(x,y,-K),t) = - 


Uoe^ 

2 

2JtP 


iCDTn iCUTp 

e • + e 


, d| dq (40) 


=’wl 


where ®wl d,s the area PABC as shown in figure 6. Since; the two- : 

dimensional solution' with supersonic leading edge' is known, the potentiajL 
can be obtained alternatively by calculating the contribution of S' ' ' ' 
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which is "the axea BCD^ the difference between the two-dimensional source 
area PAD and as shown in figure 6- Then, 


0r(x,y,+O,t) 




= - ^(x,y,+0,t) 


(^1) 


^(x,y,+0,t) = 



-iorr-. 

e 


+ e 


-icuTg 


r 


d| dTi 



- Tj) + w(|,Tl,-K),t 
r 



d| dll 


( 1 ^ 2 ) 


The variables Tg, r, and w are the same as before. In this 

development, it is found more convenient to introduce a new coordinate 
system as follows: The x'-axis passing through the fixed point (x,y) 

runs upward along the right forward Mach line. The origin is taken as 
the point where the x'-axis intersects the wing-tip side edge. The 
y ' -axis following the other Mach line runs left upward as shown in 
figure 6, The origin of the x',y' coordinate system is translated 

from the origin of the x,y system a distance where 5 can be 

shown to be 
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6 = ^ (x + Py) (I<-3) 

which, is a function of the point P(x,y) only. The transformation 
relations are 

- S ® - »y) = ^ - s) X =-|(x' + y' - S) 

y “ - 1 3 + py) = -|(l y = s(''' - y) 

-5® - I =-|(^' +V -8)>(1A) 

- I ® + ^i) >1 = |(l' - t|') 

3(x,y) 2p 

and 

^ ^ \/(-'n')(x' - I') = ^ ^Ti'd' - 3t') 

^1 = " (^5) 

""2 = ^ 



Nov and can "be expanded in terms of the nev variahles 


--UkA^- J_ 


^ f-iC0Ti')n ^ 

2_ ' r '-- = 2_A-i.r^r. 


n=0 ^ 

CO 


n=0 


n; 


-i^T 


D.=u 

vhere p = - x' is introduced. 

Similarly, 

-itDTi 




(4t) 


-IXDTg 


If attention is restricted to 1 ojTt ( < 1 and loorpl <1, e" and e ' can be 

I — I . ^ — I 

approximately represented by a few terms of the series . In the later development, five terms 
of each series ai-e used. Of course, for higher values of laJi"i| and |uiT2| more terms should 
be used. Introducing this approximation, into equation (4-2), 


CO 

o 




ia^ 




^ J J 


m 


iat 




- ^ (ii +.’i') - (t* + '■)* ■*■ ^ ’'’5^ * 


■ fu + n'r + im'ul 


«ja . 




(P + I)') .+ - ki V ^ 






(W) 
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Integrating witli respect to and Introducing a group of new 

notations, 

PS 

ic0t pp-k - ^ 

0(x',y’,-fO,t) = 


■\j2^Uaje 


dz 


JtM, 


1 / P It R 

■■ epz +:ei,z^ + e^z^ + eoz'-’ + 

ri Jo _ 

+ ggz6 + ggz8) + (t - 

(W 

„hera 7i = 72 = ^, .2 = 1, - | b = ^7l5 - | x', and 

the e*s, g’s, and h^ are given in appendix B. With ’the aid of the 
formulas in appendix C, the following integration of equation (49) can 
he carried out with respect to z : 


^(x,y, 40 ,t) = ae^ 


V ^ n=l n=0 


(|-y(x - 

*■ ' 


- 

+ ^-y(x + 

2 

Py)^ ^ oj 

n=0 

■JJ ** 


n=0 

\j-y{x ■ 

*• Py)”^ ^ ' 

n=0 

3 

■ 

■ ky)“ 

+ ^-y(x + 

%o ■ 

o 


+ Py) ^ ' 


- ky)“ 

+ \|-y3(x - 

*■ Py)3 ' 

n=l 

1 

1 7<“ 


n=0 

\[-y^(x 

+ Py)5 XI 

■rL=0 


- isy)^ 

^ + \/-y^(x 

+ Py)T' cr^o 

CD 

(m and 

n are integers) : 

in terms 


+ 


( 50 ) 


2 ind P are given in appepdix B. By definition, 

C = ^ 

^ ^t U 5x 

Differentiating equation (50) with respect to t, 

2 2 icD 


( 51 ) 


U' 


■2 at 




V‘= 


( 52 ) 



1 


Differentiating equation (50) with respect to x and adding equation (52)^ 


— / \ 1 _ iXDt 

Cp^x,y,+0;-o; = - “I Bin 

\J-y(x + py) 

3 


-1 / x + 3y 
\/x - ky 


n=l 


/ _ \n-l . ItD 

^ - “y; ■" T ^ 


1 - U02 ^(x r ky)*^"^ + ^ ^ ~ + | ^ ^3 + 

M=1 ^ ^ ^ n-Tl 


n=l 

_2 

> 

n=0 
1 




+ \/-y(x + py)3 


- ky)“‘^ + T ^ ^ 

Gr-’l il®0 


I X + fy(i + 3y)5 ^ noi^ n(x - ky)""-^ + ^ X 

n=0 n=l •' n=0 ' 


2 ^ 


n-O 


. + \/-y(x + py) 


-5,x * f ± - I -6,0 


\|-y(x + Py)^f^ 06. o) + \fy^XX7) ^ „(x - ky)“-^ + ~ + 

•'“/ ' -nol n=1 ■■'“ 


2 . 


5^ °8,n(* - V)“ 


2 ^ "S'” 


(x - ky)^ 


+ \/-y^(x + p^)3 

+ \|-y^(x + py)5 


n=l ■'■ 
2 


na0 (^ - ^X. + 

n°x - n=0 


0 + ^ 
9,1 U 


- ky)^ + i cT;^o,0 


\|XXXX f 010^0 


(53) 


LO 
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where again the 
coefficient is 


CTjjj^n'’^ are given in appendix B. The acttxal pressure 

2 S0IJI 2 

~ 'u^ ^ " u 

J_ / A M 

y2 ^t U' ^2 St U Sxy 


= Cpo(x,y,0,t) - Cp(x,y,0,t) 


(^) 


where 0 t = given in equation (4l) and CpQ(x,y,0,t) is the 

pressure coefficient of the two-dimensional case as given in reference 5- 


) 


Rectangular Wing Tip 
For the rectangular wing (k = O) equations (53) and (5^) reduce to 

4 


Cp(x,y,+0,t) = I oe^ 


sin"^ J - ^ + W-y(x + py) i>n^ + 

V ^ n^ ^ n=0 


n=0 

_1 

n=0 


2 

\/-y(x + py)3 ^ Cjj^x^ + \j-y(x + py)5 ^ d^x^ + ' 

n=0 


\|-y(x + Py)"^ ^ e^x^ + \j-y(x + By)^ f + 


n=0 

J, 

n=0 


2 

\|-y3(x + By) ^ gj^x^ + \l-y^(x + Py)3 ^ hj^x’^ + 
- ^ n=0 


\J-y3(x + Py)5 i^x^ + \j-y^(x + Py)”^ J 


(55) 


where the a's^ h's, c'b, d's, e’s, f, g's, h's^ i's, and j 
are given in appendix B. (The particular Ujjj j^'s for k = 0 are 

found in appendix B . ) 
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Average Presexore Coefficient over Rectangular Wing Tip 


For engineering purposes, an average value of over the 

entire tip influenced by the three-dimensional flow may be useful. 
(Assume unit tip chord.) 



where 23 is the wing -tip area influenced by three-dimensional flow. 


Si 


-1 - icafc 


n-x/3 


^0 


n=0 


<3y 


^0 




x+Jy 


— 

2 , 




I ^ 

1 I (jx y I cos”'^ q dq 


JO 


-1 - icct ®n 




oe 




n=0 


+ 2 


(57) 


where 


X + 23y 


q 


X 
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r - 


c/p 


dx I dy 

/o Jo 


^-y(x + py) ^ b^x^ + 
n =0 


,/-y(x + Py)3 ^ c^x^ + J-y(x + 3y)5 ^ d^x^ + 
V r,=n 'I n=0 


^-y(x + py)'^ ^ 6 j^x^ + \j-j{x + Py)^ f + ^-y3(x + Py) ^ g^^x^ + 


n =0 

n =0 

DL=0 


n =0 


\/-y3(x + py)3 ^ h^x“ + ^-y3(x + py)5 ^ + J-y3(x + Py)"^ j 

'' — ^ n =0 ^ 


( 58 a) 


Let Py = -z^. -E 3 q)£ind 1 .ng higher powers in terms of (x - 

the integral may be integrated by formula (C2) in appendix C to 


^P2 


i oe-^ dx^bj'x^ + b2'x3 + b^'x^ + bj^’x5 + b^'x^^ 


-1 - iout 




ae 


n^ ^ + 2 


( 58 b) 


The b^' s are given in appendix B- Hende, 


C = li aei^ ^ ^ c ' 

P^av 4p ^ n + 2 n 


n=0 


n=0 


(59) 


(See - appendix B. ) 


Calculations and. Graphs . - 

Owing to the complicated nature of the formulas in this section,, 
the calciilations are rather tedious. The calculations are too long to 
be included here, but the essential results are summarized in figures 7 
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to 15 . The pressiire coefficient given in this, paper is referred to a 
point on the upper siirface of the wing. Figures J to 10 give the pres- 
sure distribution on a rectangvilar wing tip of unit chord with different 

CDC 

free-stream flight Mach numbers . The frequency parameter — is taJsen 


as 0.2 in these calculations. As before the complex quantity 


JtC, 




is used and plotted in two parts - real and imaginary (e.g.^ figs. TC^’) 
and 7(t>) ) . These curves are plotted along the spanwise direction (y 
being negative when measured inward from the tip) at various percentages 


^T) 

of chord. If — along the chordwise direction is desired, some 
oe^ 

cross -plotting is necessary. It should be noted that these curves do 
not give the pressure coefficient Cp directly, but Cp may be calcu- 
lated easily from them by using equation (37a) or (37t>) or both. 


Figure 7 is calculated for a Mach number of \[2. (Note that 
figure 7 checks well with figure 5 which was calculated by the first 


method.) It is seen that 


iiC, 




is zero at the wing tip (see 


fig. 7(a)) Tor all valties of x, except at the leading edge where 


itC, 


^oe^ 


becomes 2it, which is in agreement with the two-dimensional 


R 


result along the Mach line -y = x. Taking the trailing edge (x = l.O) 


as another exanple. 


jtC, 


oe^’l 


increases from zero at the tip (y = O) to 


a value of -1.885n at -y = x = 1. Thus it is seen that for this 




case 




R 


varies veiy little along the left Mach line of the tip 


as is shown by the dotted line in figiire 7(a). The horizontal solid 
line for -y > 1.0 is the two-dimensional result'. The imaginary part 






shown in figure 7(t) is of a similar nature. Figures 8, 


'I 


and 10 show the same wing oscillating under the same conditions except 
that the free-stream Mach number is now 2, 3; and 4, respectively. A 
few points may be of interest. As the Mach number Increases, the maximum 


jtCt 


nC 


magBltuaes of UotH 
' dimensional values become more nearly constant. 


decrease, while the two- 
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In order to show the accuracy of this method, a comparison between 
the result of the present approximate method and the two-dimensional 
values along the left Mach line of the tip cone is given in* figure 11. 

The case considered is a rectangular wing tip (k = O) and the free-stream 
Mach number is teiken as ][2. The solid line shown is the approximation 

COO 

and the dotted line is the exact one. At — = 0.283, the deviation is 

ap2 

CDC 

unnotlceable . At — = O.'JO'Jj which is beyond the in^josed limit, the 

aP^ 

maximum deviation of the imaginaiy part amounts to 10 percent. The same 
curve at x = 0.8 'coiresponds to the case of x/c = 1 with lower fre- 
quency parameter = 0.707^ a't ^ = 0.5656 the maximum deviation is 

ap"^ ^ 

coc 

of the order of 2 percent. Thus beyond — — = 0.5 the approximation of 

aP^ 

the Imaginary part is not very favorable unless more terms of the series 
expansion are taken. 


There is another parameter which influences the accuracy of the 
present method, that is, the sweepback angle (a = tan“^ k) of the super- 
sonic leading edge. Figure 12 offers such a con5>arison. Here the Mach 

COC 

number is V2, — ^ = O.283, and the parameter k varies. The heavier 

aP'^ 

line gives the two-dimensional exact results along the left tip Mach 

Cp 

line, and the lighter line of the same kind gives of the present 


oe 


method . Although 


oe^^ 


does show slight deviation from the exact value 


at higher values of k, the comparison is favorable. 


In order to show some cases with high sweepback angle (although the 

leading edge is still supersonic), figure 13 gives the distribution of 

C-p . — 

— over the three-dimensional wing tip of a 45° sweepback angle at 


M = 2. For a matter of convenience in calculation, the results are 
expressed in oblique coordinates x + Py and x - y. The line 
X + Py = Constant is parallel to the left Mach line, while 
X - y = Constant is parallel to the leading edge. Here, the chord 
normal to the leading edge is chosen as unity. Therefore, x - y 
represents percentage of chord. With such a definition of the chord, 

CDC 

the chord c used in the frequency parameter — ^ should be considered 

ap2 

as this normal chord. The circles at the left lower comer of 
figure 13(a) represent the values along the left Mach line. With the 
sketch on the figure, the general features of the ciorves are clear. 
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A comparf-son of results of t his sweephack case with two- 


dimensional exact -values of 


oe^ 


•is also shown in figure l4. The 


con 5 )arison is fa-vorahle. The -two-dimensional -value of 
swepthack -wing is obtained from reference 6. 


;=„iajt 

oe 


for a 


For engineers, the pressure distribution on the three-dimensional 
wing tip is not of so much interest as the average pressiire distribution 
Cp Figure 15 calculated from equation (59) inay serve such a purpose. 

1h.e expression = -2Cp^g^^ may be considered as the lift coefficient 

for the three-dimensional wing tip- The rec-tangular "wing tip (k = O) 
at -various Mach numbers is considered. The dotted ctirves correspond to 

GiC CDC 

— = 0.2, and the solid curves, to — =0.5* The magnitudes of both 

& 3i 


the real and the ImELginary parts of . decrease with increasing 

CiC 

Mach number as mentioned before. The change of — from 0.2 to 0.5 

a 

Cp 

causes a sli^t drop in the real part of — but a greater drop in 

oeiast; 

CDC 

the imaginary part. The average pressure at — =0.2 de-viates so 

a 

little from -bhe steady case (to = O) -that it is hard to show in 
figure 15 (a) . 


The present method does gl-ve a closed expression, but it takes a 
great deal of time to calculate even a single point. 


DISCUSSION OF EESULTS 


With the above analysis, a few points should be emphasized. From 
Eward's concept of the "diaphragm,” the two approximate methods have 
been shown to give a satisfactory pressure distribution in the wing tip 
influenced, by the three-dimensional flow, if the frequency parameter is 
low and if the leading and trailing edges are si^jersonlc . 


The comparison of the pressure coefficient obtained by the second ’ 
method -with two-dimensional exact solutions of linearized theory is 


ox 


favorable, if the frequency parameter — q 1b around 0.5 and the point 
is at the trailing edge . It is expected that the first me-thod should , 
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give satisfactory results wltli the frequency parameter of the order of 
one-half. Of coiurse. If only the average pressure coefficient over the 
entire tip hounded by the left tip Mach line Is wanted, a higher fre- 
quency parameter up to the order of unity may be used. 

He effect of the three-dimensional flow on the wing tip Is essen- 
tially to reduce the pressure Increment In the area. With the approxi- 
mation given by the present theory, the pressure coefficient at the tip 
side edge Is zero except at the point on the leading edge. Of course, 
the exact theory would give finite values at the side edge, "vdilch would 
approach zero (as cs — >0) as predicted by the conical-flow theory. At 
low values of the frequency parameter, the effects of the finite edge 
pressures give negligibly small contributions to the aerodynamic behavior 
of the wing as a whole. It Is In this fact that the present methods find 
their Justification. 

The effect of Increasing Mach number Is to reduce the real and 

Imaginary parts of ^ . In magnitude. If the sweepback angle and the 

frequency parameter remain the same. If other conditions of the wing 
remain the same. Increase of the sweepback angle causes the magnitude of 

Cp 

the real part of to Increase and the Imaginary part of It to 

oe 

decrease as shown In figure 12. Increase of frequency parameter has 
the opposite effect. 

The pressure coefficient of a point Influenced by the nose of the 
sweptback wing Is given In terms of the Presnal Integral. It Is not 
difficult to carry out the numerical integration for this case. In this 
portion of the wing, no approximation is used. 

The first method of calculating Cp at the tip using the Fresnal 
Integral appears worthy of more extensive study in the future as the 
method Is very simple and useful. The second method gives closed 
expressions for Cp and Cp and it has been Investigated extensively 

In the present report. 

The present analysis applies only to cases in which the flow at the 
leading and trailing edges Is supersonic, and the nose Mach cone and the 
tip Mach cone do not intersect within the wing area. With sli^t modifi- 
cation the methods can be extended to more general wing plan forme. 

The Johns Hopkins University 

Baltimore, Md., November 15, 19^9 



ho 


MCA TN 246? 


APPENDIX A 


SYMBOLS 


a velocity of sound 

A source strength 

7-1 

b = P7t6 - — x' (equation (49)) 
/o 


coefficients of C 


p,av 


chord 


drag coefficient 


lift coefficient 


moment coefficient 


pressure coefficient 


Fresnal integral 


cos ^ (e')^ de' 


function concerning source strength 
maxi mum descending rate 


molt vector in x-dlrection 


constants expressing -wing inclination 
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M 
1 
q' 

r = 

S 

s(0 

t 

u 

u,v 
u' - 
v' 

V 

X 

x',y' 

y 


free-stream Mach number 
velocity vector 
perturbation velocity vector 


- 1)^ - (m 2 - l)[(y - -n)2 + (z - 


surface 


Fresnal Integral | / sin (e')^ (see fig. 2) 

lo 


time 

free-stream velocity- 
characteristic coordinates 
source location along u-axls 
so;irce location along v-axls 
velocity component In z-dlrectlon 
longitudinal axis 
transformed axes 
spanvlse axis 
1 / 2 . 


- (S) 

1 / 2 / v„\l/2 


Y-, = 


2ojP 


iitaM^ 


I'-e) 


z = 


2cd \ 

apjtj 


vertical axis 

i/2r 


|(v„ - v')V2 - i(u„ - u')V2 
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Z3 = 


205 


dJj / ^ 

aBn ' 


V 


L^) 


kiM^ 


1/2 


2m \V^ 
aBnj 


- 


^ ] ^V~ ^2) 


a 


a 


angle of attack or slope of streamline surface 
amplitude of angle of attack 


B = \/m2 - 1 


A 

7n 


sweepback angle 
constant 


H ) 


constant 


constant 


(b i k) 


2k 


B + k 


6 = I (x + By) 


i 

T1 

\ 

X 


source location along z-axis 
source location along y-axis 
slope of diapbragm 

amplitude of oscillation on diapbragm 


^ = i' - X ' 

I souxc.e location along x-axis (t-axis) 

Og slope of bottom wing svirface 

a,j slope of top wing surface 

Ojjj^n coefficients for Cp 



T 


time interval 





6 ) + 


r 

a 


0 


CD 

0 ) 


CD 


perturbation velocity potential 
ang\ilar velocity per second 


frequency parameter 
frequency parameter 



Subscripts: 

B bottom wing surface 

D diaphragm^ except in tbe case of drag coefficient C 

I imaginary part 

R real part 

top wing surface 


T 
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APPENDIX B 

iFOEMULAS USED IN SECOND METFOD 
Fonmilas for Use with Eqmtion (49) 


62 -= 4 - .l^icD/gC X - ky ) - aj3^72^(^ “ 


ei, = 


^31 ^2 

M 7i 


-4icu - + 2ini^-/2^(^ ~ + ■§■ 


®6 = 


Im/Z2' 
m2 W. 


_ 733/72\3 


- 20^2 + 2ioj372(x - ky) + coV 2 ^(^ “ ky)2 


I iui3 + I odV 2 (x - ky) 


’'34/^2 1 -4 

3 - 


«4 


, /y \l/2 p 

m41^) "f ^ f ^2^^ - 4y) + ^ - ky)2j 




^31 

M 


X- ky)^ 

- ky)^ 



k 6 
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where 


_ ooMc 
CD = — 5- 

ap 


2 k 

^3 = 3 -Tk = 2 k /2 


^31 " 72 


(^3 - 1 ) 


'7lf( 2 ^ -8\ 

’'32 = I^K^3 - 3 ’'3 ^ I5J 


.^2 


33 



27 ^ + — 7 
^3 5 ^^3 


_16\ 

35/ 


734 


= 7 ."^ - ^ 7 3 + 1 . 2 _ II „ 


2 ^^3 “ 3 ^3 5 ’'3 ■ 35 ^3 315/ 


- ^ 1 / 2 
’'31 = - 5 


\^/ 2 4 

7oo = I 1 I ^3 - 5 ^3 + 


32 /2 



Formulas for Use with Eqixations (50) and (53) 


■’l,! = 2 (ri 72 )^/^( 2 ) 


%2 = 2 {ri 72 )^''^li®' 2(-2 - I 731)- 
■^ 1^3 = 2(ri72)^/^;5V(-i - i 731 - g 732 - 1 ^ I) 
'’1,4 = 2 ( 7 i 72 )^/%S 53 r 2 ^(| + f 731 + §732 


• 35 1 51 ^ 1 - 

■" 192’'33"’ 6m 2 72 ■" it6 m2 72’’31 
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1 1 


1.5 


1 ^1^5 1 ’'l - 


12 ^ k ^31 i6 ^32 " 192 ^33 " 256 ^3^ 


35 21 
+ 7oo + rr^ 7-.- + 


+ ^ ^ + _I_ A. Zi - + 1 Y ^i\^ 

12M2 72 ^ M^ 72 192 7 2 ^^2 l60M^\72/ 

^^2,0 = 2(7i)^/2(-2) 

'^2,1 = 2(7i)^/^Jiii72(2 + | 731) 

^^2,2 = 2(71)^252^2^^1 + i 731 + I 732 - ^ 


2,3 



,T3-y_3/_l _ 3 


+J^Zl+5_iZl- 

3 ¥^31 8^32 192^33 6m^^2 ^m2 72^ 

1 1 5 35 21 ^ 

■I 2 " ¥ ’'31 " ^32 - ’'33 ■ ^ ’' 31 ^ 

1 /7i\2 



^2,l^ = 2(7i)^/^ajV2^ 

^Il + ^Zl^ +JLl_Zly + 

12M^ 72 72 £92 7g 

^3,0 ^ ^(’'1) ^^^^2(731) 

"3,1 = 2(7i)V2iD272^(r3i + ^ 732) 

'3,2 ' ^31 - s r32 - ^ ’■33) 

^3,3 “ 2(7 i) 1/%V2‘‘(-| 731 - i ^32 - ^ ’■33 - ife ’’31.) 

74.0 = 2(7i)^/^^72^(j 732) 

74.1 = 2(7)^''^is3r23(-i 732, -. 733) 

74.2 ' 2(7)^/V72‘^(-| 732 - 4 733 - rfe ’■34) 


“5,0 = 

“5,1 ■ ’'33 - k hh) 

“ 6,0 = 2 (?)^/^a>^r 2 '*(-^ r^C) 


“7,1 = 2(71)3/2 i £62(1 


°7,2 - H^i)' 


3/2 Z2 ^ _ 




3 2k '31 


) 


'^7,3 = 2 (^ 1 ) 


3/2 Z2_ 5,4(_i _ 2-y^^ _ J.y _ i. i_ Zl 

80 72 




' 6 ~ 2k ^31 ~ 96 ^32 


. 3/2 iisCitl 


“ 8,0 = 2(71)- ^2 -, 3 ^ 


“8,1 = 2(71)3/2 731)1553 

“B,. = 2(7i)3/2/j6f(.|.^73i-^73,-|,.^|) 

“ 9,0 = 2 ( 71 ) 3/2 ^ 1 S 3 (.| i^^) 


“ 9,1 = 2(71) 


3/2 X mH-i y 


U‘ 


7 ^ 


1 1 


3 ''31 60 ^32 50 72 


^ 10,0 ^(^ 1 ) 


3/2 l2 ^k/2 


M“= 


^ Jl 1 Zi 

15 32 75 72 


a - 2(y m 2 Z£/ i Zi^ 

^2,1 - 2(71) j^2 ("6 72 ] 
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^ 5 - 

" W ^ ^31 


2,3 


M' 


1 ’'1 ^ 5 n 


7 ^3 


FormulaB for Use. with Equation ( 55 ) 
iCD 

^1 " ~u %1 ■*■ ^%2 

1X13 - 

^ - Tf %2 + 301^3 

^3 ~ IT %3 ^ 1 , 1 ^ 

1X13 _ 

H = u 

1X13 

^5 = ir %5 

, 1x13 3 

^0 = '^ 2,1 + - jT %0 + 2 *" 3,0 

’=1 = =<’ 2,2 ■>■ f “ 2,1 + i %1 + 

'’2 “ 3 c > 2,3 + ^ » 2,2 + I + 5 ^ 

, 1 lil3 3 

= ‘*'' 2,11 + T ■’ 2,3 * i ” 3,3 * “ 2,3 

, i03 

= T, 

1X13 *5 

^0 - IT ^ 3,0 2 \o 


1 

160 
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ioi 5 

= 203^2 IT ^^ 3,1 ■*■ 2 \i 


_ iCD 5 

'"a = 303^3 + cT 3^2 + 2 %2 


ioj 

'"a “ u" %3 


J , icD . 7 „ 

<^0 = %,i ^ If ^^ 4,0 2 ‘^5,0 


iCD 


^1 - ^%2 IT %1 2 %1 


^2 - If ^^4,2 


itD 9 

®0 = ^ If %0 ^ 2 %0 

im 


ei — 


U "5,1 


^ _ ICD - 

^ u ^6,0 


«o ~ ""7,1 2 %o 


iOD 


Si - 2 c^ 7,2 U ^ 7,1 2 %1 


icD 3 

g2 = 30y^3 + -jj- Oj^2 2 %2 

iXD - 

®3 “ U ^^7,3 

, io> 5 

^ If %0 2 ^^ 9,0 

, iuj „ 5 

^1 “ IT '^S,! ^%2 "^2 °^9,1 


X, icJO 
^ “ U %2 
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^0 ^ 9,1 U ° 9,0 ;*■ 2 ‘^ 10,0 


'‘■I - U "^9,1 


J = O' 


10,0 U 


Formulas for Use with Equation ( 55) 5 cd = ^ = — 

P afi 


%1 - p 
^1,2 = 


_ 1 =2/ 2 1 
^^1,3 - P \“3 " '3m2 


1^4 ■ p 


1=4/1 , 1 ^ 
%5 - P "'{so •" iom2 

<' 2,0 = 

= P"^/2 i^2) 

0-1/2 =2/2 1 

CTo 2 “ ^ 0> ( o r 

^ ,\3 3M" 



,^3 = rV2.a3(.i,^ 


»2,4 = P-"/" 


%0 
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®3,1 = ^(-|) 

“ 3,2 = ^{j) 


3,3 


= 6-1/2 


%,0 = 

%1 ■ ^(-5) 

“It, 2 “ 

“5,0 = iS3(j§^) 

^5,1 = 

“6,0 = “‘‘(w) 


7,2 


. el/2 4. jgS/'.i 


M^= 


'(-I) 


'7,3 = „2 “X-y - 


„ „l/2 1 =2/8 

%o = P ^“(9 

= bV 2 i_ ,= 3/2 


^8,1 = P 




_ p.1/2 1 =k( 4 1 

8,2 P „2 “ "15 


\ 15 3 om 2 



NACA TN 2467 


53 




■^9,0 
cT9 




15 y 


-1/2 1 =l^/l6 1 \ 

^ v5 25M^/ 


10,0 


= pV2 4§^/'^^ 8 




i525 


75M^ 




= P 


- 1/2 k _ = 2/1 


M‘ 


.2 “ 13 


^(1) 


02_^3 - 3 


= G " V 2 JL m ^/ j ^ 


80 M^, 


Formulas for Use with Equations ( 58 ) and (59) 


to' =0 


bi' = w2 


p U 


, _ 1 cd/2 =\ 

! p U\3 “) 


zi ^2/1 
P U 


03“ Hr + 


12M" 


, -1 CD = 3/1 

' = ^ (D-^f 


p U \^30 20W 


S = H(®'‘)(iio * * 



1 . 


'o 232 
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'_1 + 1 l! _ -A_\ 

^ 8 6 i6m2 / 


i. _ X P£ + _ 1 . XV 

30 2k ' 20M^ 


_1 1_ V ^ _1_ _ J:_ J|£ + 1 \ 

llA 120 l*Bl^ 80 isitf?-/ 


' 1 i£ ^_1^ ^ 

720 if. 2ko i^ 1920 m8 j 



APPEHDIX C 


INTEGRALS USED IN SECOND METHOD 


fi-X 


(a 2 . ^ 2 ) 1/2 


/ 2 3 

^ ^ ^ (2n - 2r + 2)!n!(n - l) I 


( Qn - - T* 4- 1^?fn - • 0( y'—'\\ Oti_C>t'4-T 

^ a^'* + 


( 'yr, _ 1 \ I 


2(n!)(n - 1)1 


Ol-l "1 V 

a Bin — 
a 


(n 1) 


(Cl) 


- x2)1/2 ^ _/2 _ x2)3/2^ 1 (2n - l) I (n - r-H) 1 (n - r + l) 1 ^2(r~l)^2n-2r+l 

1{.^“^ (2n-2r + 2)!(n+l)J(n-l)! 


1 (2n-l)!a^ F / ? p^l/2 p , _i x 

lt.n (u + i)l(n- 1)1 ^ a 


(n^l) (C2) 


1 ^ “02 (^^ + = 2ti^^^(|i + n) - I ti3/2 


(C3) 


ui 

Ul 


Nl VCJVH 
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+ t))^ - I ti3/2(^ + t^) + ^ (c4) 


dT] 


ri)3 = 


2tiV2(p. + t|)3 _ 4tj3/2(^ + tj)2 + ^ t]5/2(ii + t^) - 


35 ^ 


(C5) 


dTi 

;i72 


(\i + Ti)^ = 2ti^/^(ii + n)^ - ^ + t))3 + ^ + Ti)^ - 




(c6) 


dn + n) = I + n) - ^ 


(C7) 


dTl Tl^/^(n + Ti) 2 = I Tl3/2(^i + t])2 _ A ti5/2(^ + ^ ^7/2 


0 
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PigiiT^ 1-— Division into soncs of stresinXiiis suri^ce of & >rLng iuS 

tip diaphragms in supersonic flov. 


i9t(S Ml VDVM 




ON 

o 





X 

Figure 3 .- Wing areas and and diaphragm area Sjj which lie 

in forward Mach cone of a point (i,y) on wing surface. Two coordinate 
aysteme (x,y) and (u,v) shown. 
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rigure 4.- Wing areas W„i and W^p -which lie in forward Mach cone of a 
point (x,y) in nose section of a wing -with aupersgnlc leading edge. 

Two coordinate systems (x,y) and (u,t) shown. 


On 


KACA Tn 21*67 



62 


NACA TN 21*67 



Figure 5*- Pressure distribution along trailing edge of a rectangular 
wing of unit cboid with ^ = 0.2, M = ][ 2 , and x = 1. Calculated 
with first method (Fresnal integral). 



>1) 


x' 





the potential at P(x,y) 
.osen for calculation in 


ON 
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(a) Real pairb of 


Figure 7 • - Real and imaginary parts of 


along -y at various 


lues of X for a rectangular wing tip (k = O) of unit chord with 
= 0.2 and M = ^2 (used for calculating pressure distribution 
on wing tip influenced by tbree-dimensional flow) . Second method. 











(a) Eeal part of 


jtCp 

Figure 8.- Heal and imaginary parts of — — — along -y at various 

ae 

values of x for a rectangular wing tip (k = O) of unit chord with 
= 0.2 and M = 2. 
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0 .2 .4 .6 .8 JO 

X ^ percentage of chord 



("b) Imaginary part of 




ae 


lent 


Figure 11.- Con5>arlson of results of second approximation wltli two- 

dimensional solutions at left tip Macli line of rlgjit rectangular wing 

tip (k = 0) along chord at = 0.283 and ^ = 0.707 with 

ap^ ap^ 
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(a) Real part of 


ae 


lost 


Figure 13.- Distribution of 




cue I 

values of x + Py for a 45° sweptbaeik -wing tip at M = 2, — = 0.4, 

and k = 1. 


icut 


along constant x - y at various 
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("b) Imaginary part of 


ae 


lOJt’ 


Figure 13 •- Concluded. 











(a) Real part of 


Figure 15- - Average pressure distribution over entire rectangular iTlng 
tip (k = 0) influenced by three-dluenslonal flow as it varies with 

Mach number for ^ = 0.2, and = 0.5* 






